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Abstract. The excitation mechanism of low degree 
acoustic modes is investigated through the analysis of the 
stochastic time variations of their energy. The correlation 
between the energies of two different modes is interpreted 
as the signature of the occurrence rate of their excitation 
source. The different correlations determined by Foglizzo 
et al. (1998) constrain the physical properties of an hy¬ 
pothetical source of excitation, which would act in addi¬ 
tion to the classical excitation by the turbulent convection. 
Particular attention is drawn to the effect of coronal mass 
ejections. The variability of their occurrence rate with the 
solar cycle could account for the variation of the correla¬ 
tion between IPHIR and GOLF data. Such an interpreta¬ 
tion would suggest that the mean correlation between low 
degree acoustic modes is at least 0.2% at solar minimum. 

Key words: Sun: oscillations; activity ; flares; magnetic 
fields - methods: analytical, statistical 


1. Introduction 

The theoretical mechanism of excitation of solar acous¬ 
tic modes by the turbulent convection is well documented 
(Goldreich & Keeley 1977; Goldreich & Kumar 1988; Ku¬ 
mar & Goldreich 1989; Balmforth 1992a,b,c; Goldreich 
et al. 1994). High frequency oscillations were interpreted 
in different ways (Kumar & Lu 1991; Goode et al. 1992; 
Restaino et al. 1993; Vorontsov et al. 1998), in order to 
address the question of the depth of the excitation source. 
Local observations led Rimmele et al. (1995) and Espagnet 
et al. (1996) to identify the excitation of 5-minute oscil¬ 
lations with acoustic events occurring in the downflowing 
intergranular regions, rather than overshooting granules. 
Considering granules of 1000 km diameter renewed on a 
time scale comparable with their turnover time (8 min), 
there are about 5 x 10® excitations by solar granules per 
damping time (3.7 days). Even if the number of efficient 
excitations is smaller by a factor 100 (Brown 1991), each 
low degree mode is stochastically excited so many times 
per damping time, that modes with different frequencies 
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are expected to have uncorrelated energies. In particular, 
the fitting procedures used to determine the frequency, 
linewidth, and splitting of p modes often rely on the sta¬ 
tistical independence of neighbouring modes (e.g. Appour- 
chaux et al. 1998). 

Baudin et al. (1996), however, noticed a possible corre¬ 
lation between low degree p modes in IPHIR data. By 
measuring the mean correlation of p modes in IPHIR and 
GOLF data, and determining their statistical significance, 
Foglizzo et al. (1998) (hereafter F98) suggested the ex¬ 
istence of an additional source of excitation during the 
IPHIR period (end of 1988). 

The goal of this study is to use the observed correlations as 
a constraint on the physical properties of this hypotheti¬ 
cal additional source. The theoretical relationship between 
the occurrence rate of this source, its contribution to the 
energy of low degree p modes, and their correlation co¬ 
efficient, is established in Sect. 2. Available observational 
constraints are reviewed in Sect. 3. The possible role of 
comets. X-ray flares and coronal mass ejections (CMEs) 
is discussed in Sect. 4. 

2. Correlation coefficient between two modes 
excited by a single mechanism 

2.1. The energy of a mode seen as a random walk 

Let us consider impulsive excitations, distributed over the 
solar surface. The exciting events, indexed by fc, occur 
at the random poissonian time tp. They correspond to 
a radial velocity perturbation Vk(j) localized in a cone 
of angle ak around the random direction {0k,4>k), with a 
radial extension Srk. Vk{r) is described by an amplitude 
Vk and dimensionless shape functions gk and hk as follows: 

Vfe(r) = gk{a)hk{r)vk, (1) 

hk{r) ~ 1 for > r > R.^- 5rk, (2) 

gkia) -- 1 for a < ak. (3) 

a(0, Ok^ff — cjik) > 0 is the angle made by the direction 0, cf 
with the direction Ok,(j)k- Let Mk be the mass of gas in the 
volume defined by the functions gk and hk , so that the ki¬ 
netic energy of the perturbation is £k = MkVk/‘^- The 








2 


T. Foglizzo: Are solar granules the only source of acoustic oscillations ? 


mean interpulse time is defined as At =< tk+i — tfc >■ 
We assume that each impulsive event triggers free oscilla¬ 
tions of the set of p modes. Each perturbation is projected 
onto the basis of eigenvectors. The damping time de¬ 
pends in principle on the mode considered, but is nearly 
constant (~ 3.7 days) for the modes considered by F98, in 
the plateau region between 2.5mHz and 3.5mHz (linewidth 
of 1/rHz). 

In appendix A it is shown that the total energy of the 
oscillations associated to the frequency tUni + mfl is: 


Klmit) = 


^nlm 

nim 


tk<t 

9Li0k)h^uiMjvk, 

{ujni + mVL)tk + m4>k. 


( 4 ) 

( 5 ) 

( 6 ) 


where the real function = gf_^ depends on the an¬ 
gular shape of the excitation projected onto the spher¬ 
ical harmonics, and depends on the radial shape of 
the excitation projected onto the radial part of the 

eigenfunction: 


gLi^k) = qir. 


dcf) / d(cos0)P;'™'(cos0) 


X cos{m(j))gk{0,9k,(l)), 


= 


Mi 


' Rt-Sric 


pr <i(r)/ife(r)dr. 


( 7 ) 

( 8 ) 


The relationship between the spherical harmonics and the 
Legendre associated f unct io ns P ™ through the constant 
qim is recalled in Eqs. (^)-(^). The exponential damp¬ 
ing in Eq. (i can be schematized as a term selecting the 
finite set of excitations which occurred within one damp¬ 
ing time before the time t: 




2^ flnZme 
0<t-tk<Td 




(9) 


The series of phases tOnitk and (pk are independent random 
variables uniformly distributed in the interval [0, 27r]. The 
series defined by Eq. is therefore also uniformly 

distributed in [0, 27r]. Ignoring the academic case where 
the ratio of the frequencies is a simple rational number, 
^nim ^^*4 ^n'Vm' 4)e Considered independent for two 
different realistic modes. 

As a consequence, the energy of the wave is interpreted as 
the squared length of a random walk in the complex plane. 
Each step of this random walk is defined by an amplitude 
^ phase (negative values of the amplitude 
can be converted into an increment of the phase). The 
number of steps N is the number of excitations, over the 
whole surface of the sun, within one damping time: 



( 10 ) 


The central limit theorem ensures that the real and imag¬ 
inary parts of the complex sum (|^) converge towards in¬ 
dependent normal distributions if N is large enough, re¬ 
sulting in an exponential distribution of energy (see also 
Kumar et al. 1988). 

According to this random walk interpretation, two dif¬ 
ferent modes excited by the same series of events must 
have correlated energies. This correlation should approach 
100% if the number of steps of this random walk is small, 
i.e. if the interpulse mean time At is longer than the damp¬ 
ing time Td- 

The independence of the phases and how¬ 

ever, makes the correlation decrease to zero when the num¬ 
ber of excitations increases. 

In particular, the random walks associated with the two 
components Uni ± rnO of a mode n, I have the same series 
of amplitudes, but have two independent series of phases. 
This leads us to expect that waves travelling in opposite 
directions do not have the same instantaneous total en¬ 
ergy (except on average), although they are excited by 
the same events. 


2.2. Theoretical correlation between the energies of two 
modes 


The mean energy < Enim > of the mode nlm is directly 
proportional to the mean energy < enim > received from 
each excitation (appendix B): 


^ Enim ^— 2 ^ ^nlm ^ ■ ( 41 ) 

The correlation between two modes nlm, n'l'm' 

excited by the same source described by Eq. (||) is derived 
in appendix B: 


'Tfi'l' rr 
■'nlm 


^ ^nlm^n'l'm' ^ 


< e , > 2 < e 

nlm n 


>2 


X [l + N 


^ O-nlm ^ 


l + N 




< 


> 


( 12 ) 


2.3. Latitudinal distribution and sizes of the excitations 

The ratio < > / < &nim >^, related to the spread of 

the distribution, is called kurtosis (it is sometimes defined 
with an additional constant —3 which we omit here). 

< ^ o^ 

< enim >2 < gfn^hhS >2 ' 

The kurtosis of Cnim is partly produced by the kurtosis 
of £k (independent of nlm), but also by the projection 
of the spatial distribution of excitations {0k, cxk, Srk) onto 
the eigenfunctions nlm. 

The random variations of the size ak,Srk of the excita¬ 
tion play a negligible role as long as it is smaller than the 
wavelength of the mode nlm. For the modes I = 0 and 
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I = 1 analysed by F98, we restrict ourselves to perturba¬ 
tions such that Ofc < TT, and Srk is smaller than the depth 
of the first radial node of the eigenfunction Vni, and thus 
neglect the random variations of ak,Srk. 

Let us estimate the kurtosis of the distribution gimi^k), 
due to the projection of the latitudinal distribution of ex¬ 
citations on the mode nlm. From Eq. (|^), the function 
gim{&k) is independent of 6^ for the radial mode I = 0: 


^ 9o,o > 
< 9lo >" 


(14) 


By contrast, the mode Z = 1, m = ±1 is more excited 
by equatorial events than polar ones. For small scale ex¬ 
citations (ofc = 0) distributed uniformly over the sphere, 
Eq. (Q) implies: 


< gt.±i > ^ 6 
<5?.±i>" 5' 


(15) 


The kurtosis of gf should reach a value even closer to 
unity if the excitations are distributed in an equatorial 
region, like CMEs at solar minimum or big flares. 

Let us approximate the kurtosis of enim as the product of 
the kurtosis of £k by the kurtosis of gf^{dk)- 


< ^ <9L > 

< enim < £ >2 < gf^ >2 ’ 


(16) 


Using this approximation in Eq. (p^, the correlation be¬ 
tween low degree modes I = 0,m = 0 and I = l,m = ±1 
excited in the random direction {9k,<pk) by a small scale 
excitation (a^ = 0) is: 


^n',0,0 

^n,0,0 


^n,0,0 




(l + iVeff)-\ 

(^) + 



5NeS 

6 


-1 



(17) 

(18) 

(19) 


where we have incorporated the kurtosis of £k into the 
definition of the effective number Ngg of excitations per 
damping time: 


( 20 ) 

N 

A'eff = (21) 

The correlation displayed in Fig. appears to be hardly 
sensitive to the geometrical content of the function 
since all curves are similar within 10%. Error bars of 
the observationally measured correlations being typically 
larger than 0.1, Eq. ( 0 ) shall be considered accurate 
enough for the purpose of this study: 


^n' ,l' ,m' 


il + N,g)-\ 


( 22 ) 


This is equivalent to approximating the kurtosis of enim 
by the kurtosis of £. 



Fig. 1. Theoretical correlation between two modes I = 0 
(full line), I = 0 and I = l,m = 1 (dotted line), and 
I = l,m = ±1 (dashed line), as a function of the number 
of effective excitations per damping time (ck^ = 0). Hor¬ 
izontal lines indicate the measured correlation in IPHIR 
data and the upper bound found in GOLF data. 


2.4- Case of two excitation mechanisms superimposed 


F98 used a one parameter model of two sources of excita¬ 
tion, such that each mode contains a fraction A of energy 
common to all the modes. They assumed that both sources 
produce exponential distributions of energy. This simplifi¬ 
cation led to the simple relation A = This simplifica¬ 
tion, however, does not apply to the case of a correlation 
produced by rare impulsive events. If iV <C 1, the energy is 
damped to zero except during isolated pulses of duration 
^ Td/2. The distribution E thus follows a Bernoulli statis¬ 
tics, with Var(£’)/ < E 1 (Eq. B^) instead of 1 for 
an exponential distribution. The occurrence rate of the ex¬ 
citation mechanism therefore appears to be a key param¬ 
eter which cannot be neglected. Let us consider two exci¬ 
tations mechanisms acting simultaneously on two modes 
n, I, m and n', V, m'\ 

(i) excitation by granules, with a high occurrence rate 

{Ni » 1). 

(ii) excitation by additional events with a total acous¬ 
tic energy occurring on average fVad times per damp¬ 
ing time, and contributing to a fraction \nim of the power 
of the mode nlm. 

We show in appendix C (Eq. p.7D that the correlation 
between the energies of two oscillators excited by these 
sources is the same as in Eq. but multiplying the 

kurtosis of Cnim by A^;^. The latitudinal distribution and 
sizes of the exci tations play a relatively small role accord¬ 
ing to Sect. The effective number iVeff of excitations 
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per damping time depends on the kurtosis /3ad of fad: 


Pad 


(23) 


Making the additional assumption that the fraction 
AraZm = A varies little among the modes considered by 
F98. The mean correlation between these modes is: 


C 


1 + ^) 
A2 j 


-1 


(24) 


A significant correlation can thus be produced by a source 
representing only a small fraction A of the total power of 
each mode if Ngs is small enough. 

Interpreting the observed correlation C as a consequence 
of an additional source of excitation therefore requires it 
to contribute to a fraction A of the total power, deduced 
from Eq. (|^: 


A 


{^r 


(25) 


According to Eq. the fraction A of power coming 

from the additional source is related to the mean acoustic 
energy input < Cad > into the mode nlm: 


X _ A^ad ^ Cad ^ 

“^ <E> ' 

Eq. (^ ) becomes: 

^ Cad P* — 


C <E> 


/3ad 1 - C A 




1-C 


< E > 


/V2, 


(26) 


(27) 

(28) 


Let us assume that the occurrence rate A^ad varies from 
A^min to A7max with the solar cycle, while the properties 
of the exciting events (< Cad >,/3ad) remain unchanged. 
We use Eq. (|^) to express the relationship between the 
minimum and maximum correlations CminjCmax : 


In appendix C, Eq. ( p^ is rewritten as the fraction of the 
mean acoustic energy < fad > which must be injected into 
each mode nlm in order to produce the observed correla¬ 
tion: 


^ Cad ^ 
^ fad ^ 


/ad (//;fV(f)df)' 


1-C 


< E > 


.(32) 


This fraction is therefore minimal when the range of ef¬ 
ficient excitations fi, £2 contains the range of energies 
where the product £^p{£'-£) is maximal. Using in Eq. (28) 
the occurrence rate A^t and kurtosis /9 t of the distri¬ 
bution of total energy fx instead of the distribution of 
acoustic energy fad leads to a lower bound of the ratio 
^ Cad ^ ^ fad ^ ■ 


< Cad > ^ 1 2 f C I ^ <E> 

< fad > /ad \ ^ ~ ^ \ < St > 


(33) 


For a given excitation mechanism one can estimate the 
number A/” of p modes with a wavelength longer than the 
size of each exciting event, which receive a comparable 
amount of energy from the source. Considering that the 
low degree modes in F98 belong to this set. 


^ fad ^ ^ ^ Gnlm A/” < Cad ^ ■ (^4) 

nlm 

From Eqs. (|^)-(|3^) we deduce a constraint on observable 
quantities, which shall be useful in Sect. ^ in order to 
discriminate between possible excitation mechanisms: 


2 f C <E> ^ 1 

A^|/3|U-C/ <£t> ~ Af Af 


(35) 


3. Observational constraints 


A^max 1 fmax / ^ E > 

A^min f max \ ^ E ^ max 

2.5. Ejficiency of the generation of acoustic waves 

Since the distribution of acoustic energy fad produced by 
the additional source is not directly observable, we are 
led to assume that it resembles the observed distribution 
of total energy fx of the additional source. Let px(f) be 
the density of probability of the source of energy fmin < 
£t < fmax, occurring A^x times per damping time. Let us 
assume that the excitation of low degree modes is efficient 
only in the range fi < fx < f 2 ) inside which the efficiency 
/ad = fad/fx is constant: 

fad = /adfx for fl < fx < f2, (30) 

fad = 0 otherwise. (31) 


3.1. Exponential distribution of individual p modes 
energy 

The distribution of energy of low degree p modes agrees 
reasonably well with an exponential distribution (Chap¬ 
lin et al. 1995, 1997 for BiSON data, F98 for IPHIR and 
GOLF data). Chaplin et al. (1995, 1997), however, no¬ 
ticed significant deviations in the high energy tail of the 
distribution, which could be due to an additional exci¬ 
tation mechanism. They estimated that the probability 
that such deviations occur by chance during the period 
of observation is only 0.1%. Among 22512 events covering 
18331 hours of BiSON data from 1987 to 1994, Chaplin 
et al. (1997) found 51 events above 6.5 times the mean 
energy, whereas less than 41 would be expected in 90% of 
the cases for an exponential distribution. Crudely speak¬ 
ing, about 10 events are unexpected in the distribution, 
suggesting A^ad > 0.05. 
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3.2. Solar cycle variations of the total power of low 
degree p modes 

The typical energy of a low degree p mode in the range 
2.7mHz < V < 3.4mHz considered by F98 is < if >~ 
8 X 10^”^ ergs (Chaplin et al. 1998). The velocity damping 
time being ~ 3.7 days, the flux of energy required to 
excite this p mode is2<if> /tj, ^ 5 x 10^^ ergs s“^. 
According to Libbrecht et al. (1986), the total energy in 
all the p modes (about 10^) is 10^"^ ergs within a factor 
10 . 

Anguera Cuban et al. (1992) and Elsworth et al. (1993) 
measured a global 30% decrease of the power of low de¬ 
gree p modes at solar maximum. This decrease seems to 
preclude a high value of the fraction A of the total power 
which some additional source could contribute to. Nev¬ 
ertheless, the physical mechanisms by which the p mode 
power might decrease, such as damping by active regions, 
or modification of the properties of the convection, have 
not yet been quantitatively estimated. This might be ef¬ 
ficient enough to dominate the energy input due to an 
additional source. Moreover, the measurement of the am¬ 
plitude of global p modes is influenced by the presence 
of active regions covering a significant fraction of the so¬ 
lar surface at solar maximum (Cacciani & Moretti, 1997). 
Given these uncertainties, no hrm constraint can be de¬ 
duced from these observations. We shall consider ’’likely” 
a fraction A < 30%. 


3.3. Observed correlations 

F98 determined the mean correlations between the en¬ 
ergies of low degree p modes at two different epochs: 160 
days in 1988 near solar maximum using IPHIR data {I = 0, 
19 < n < 23 and Z = 1, 18 < n < 23), and 310 days in 
1996-97 near solar minimum using GOLF data {1 = 0 and 
/ = 1, 17 < n < 25). The mean correlation coefficient C 
they measured is: 

(i) C = 10.7 ± 5.9% in 1988 (IPHIR data), 

(ii) C < 0.6% in 1996-97 (GOLF data). 

According to F98, the probability that the correlation 
measured from IPHIR data could occur by chance is 0.7% 
if the modes were independent. F98 also rejected the pos¬ 
sibility of an instrumental artefact by checking that the 
noise of IPHIR data at different frequencies is not corre¬ 
lated. 

If the granules were the only source of p-mode excita¬ 
tion, the correlation would be less than 10“®% according 
to Eq. (p^, i.e. well below the detection limit. The actual 
limit of detection, of order 0.6% (F98), corresponds to at 
least one exciting event every 29 minutes, on average. If 
due to a single mechanism of excitation, the correlation 
measured with IPHIR data would correspond to an excit¬ 
ing event every 12 ± 7.2 hours. 

Fig. 1^ shows the relationship between A^eff and A required 
by the correlations observed in IPHIR data, and the upper 



Fig. 2. Relationship between the effective number A"eff of 
excitations per damping time and the fraction A of the 
mode power due to the secondary mechanism of excita¬ 
tion, in order to match the observed correlation in IPHIR 
and GOLF periods. The value of the correlation is indi¬ 
cated on each curve. 


bound set by GOLF data. 

The correlation measured by F98 in 160 days of IPHIR 
data imposes that A’ad > Td/I60 = 0.02 in the IPHIR 
period, which is coherent with the fraction of abnormal 
events in BiSON data. 

Applying Eq. (^)-(^) to a typical mode of energy 8 x 10^"^ 
ergs and damping time Td = 3.7 days, the constraint 
C > 4.8% leads us to look for a mechanism occurring at 
most a few times per day in 1988, with an acoustic energy 
input of at least 10^® ergs per mode: 


< 7.1 X 1 


> 6.7 X 1 



0.3 

0.3 \ nc 
— I X 

A 


(36) 

(37) 


4. Some candidates: comet impacts, X-ray flares 
and coronal mass ejections 

4-.1. Comet impacts 

Sungrazing comets of the Kreutz group are small rem¬ 
nants of an earlier passage of a larger progenitor comet. 
Among the 10 sungrazers observed by SMM, 4 appeared 
during the 160 days IPHIR campaign (MacQueen & St. 
Gyr 1991). A total of 59 sungrazers were observed by 
LASGO aboard SOHO. Even with an occurrence rate of 
0.1 per day, C > 4.8% in Eq. ( p^ ) would require an en¬ 
ergy input of X 10^^ ergs per event and per mode. 

Given the small scale of the comet, all 10^ p modes should 
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Table 1. Energy, momentum and occurrence rate of a free 
falling comet, a X-ray flare, and a CME (Zarro et al. 1988; 
Crosby et al. 1993; Hundhausen 1997) 



comet 

flare 

CME 

kinetic energy (ergs) 

average 

- 

~ 10®° 

6.7 X 10®° 

maximum 

2 X 10®^ 

~ 10®® 

4 X 10®° 

kurtosis P 

- 

- 33 

4-20 

momentum (g cm s ^) 

average 

- 

~ 10°^^ 

~ 10°® 

maximum 

6 X 10’^® 

- 

~10°" 

occurrence rate (day ^) 

solar minimum 

< 0.1 

1 

0.2 

solar maximum 

< 0.1 

19 

3 


receive the same energy input as low degree modes, which 
corresponds to a total acoustic energy of ~ 10 ^“^ ergs per 
event. This already exceeds by two orders of magnitude 
the kinetic energy of a comet with the mass as Halley’s 
hitting the solar surface; not to mention the variability 
issue between IPHIR and SOHO periods, nor the difficult 
question of the efficiency / of the energy transfer (ad¬ 
dressed by Gough 1994; Kosovichev & Zharkova 1995). 


4-2. X-ray flares 

4.2.1. Observations of waves generated by flares 

Although sunspots are known to absorb high degree acous¬ 
tic waves (Braun et al. 1987, 1988), observations of the 
waves generated by a large solar flare on 24th April 1984 
by Haber et al. (1988a,b) showed a 19% increase of power 
in outward travelling wave, dominating the sunspot ab¬ 
sorption. Nevertheless, Braun & Duvall (1990) also ob¬ 
served the wave emission from a flare on 10th March 1989, 
and found an upper bound of 10 % power increase, if any. 
Very recently, Kosovichev & Zharkova (1998) analysed the 
shock wave produced by the ” fairly moderate” flare of 9th 
July 1996, observed by MDI aboard SOHO. The wave am¬ 
plitude is associated to a momentum a factor 30 smaller 
than the one expected from their theoretical model. This 
unexpectedly high efficiency led them to conclude that the 
seismic flare source might be located in subsurface layers. 
The effect of flares on low degree modes is however less 
clear. The high energy events of BiSON data, appearing 
above the exponential energy distribution of low degree 
modes, do not seem to be correlated neither with the 
sunspot number, nor with the strength of X-ray flares 
(Chaplin et al. 1995). Using the sum of normalized en¬ 
ergies of low degree modes in IPHIR data, Gavryusev & 
Gavryuseva (1997) found an anticorrelation between big 
pulses in p modes and the mean solar magnetic field but 
no correlation with the sunspot number. 


4.2.2. Theoretical estimates of p mode excitation by 
flares 

The theoretical estimate of the energy transmitted from a 
flare to p modes can follow three different approaches: 

(i) the region surrounding the flare is heated and ex¬ 
pands vertically on a time scale short enough to communi¬ 
cate momentum to the atmosphere below it (Wolff 1972). 
This process could extract 10^® ergs of acoustic energy 
from a 10®^ ergs flare (/ ~ 10“^). According to Wolff 
(1972), the error bar in this estimate is at least a factor 10. 

(ii) the plasma flows down towards the foot points of 
the field lines, and hits the solar surface, thus communi¬ 
cating momentum to it. This approach was followed by 
Kosovichev & Zharkova (1995) who found a smaller effect 
than that observed by Haber et al. (1988a). 

(iii) the pressure perturbation associated with the re¬ 
structuring of the magnetic field in the flare region might 
be more effective, as suggested by Kosovichev & Zharkova 
(1995). 


4.2.3. Energy distribution of flares 

The total energy released by a flare (less than 10^^ ergs to 
1033 ergs) is usually computed assuming that the observed 
hard X-rays are produced by bremsstrahlung from a dis¬ 
tribution of accelerated electrons impinging upon a thick 
target plasma. Using the Hard X-Ray Burst Spectrome¬ 
ter on the Solar Maximum Mission satellite from 1980 to 
1989, Grosby et al. (1993) showed that the flaring rate 
varies by about a factor 20 between the solar maximum 
and minimum, while the energy distribution remains con¬ 
stant. Applying Eq. (^) to the correlation observed by 
IPHIR, with the 30% variation of the p mode total power, 
suggests the following correlation Cmin at solar minimum: 

0.15% < Cmin < 0.59%, (38) 


which is compatible with the upper bound obtained from 
GOLF data. 

The occurrence rate of flares follows a power law distri¬ 
bution p{£) ~ £“ 3 " against the total flare energy with a 
slope 7 < 2 , indicating that the largest flux of energy oc¬ 
curs in rare big energy events. Let Cmin ^ Cmax be the 
range of energies over which the power law distribution is 
observed, its mean energy < Cp > and kurtosis /Jp are: 


< Cp > 

Pf 


7-1 

2-7Uminy 

(2-7)^ /gmaxV 
(3-7)(7- 1) V'S^min/ 


(39) 

(40) 


According to Crosby et al. (1993), about 13 flares occurred 
per day in 1980-1982 (solar maximum), with a slope of the 
energy distribution 7 ^ 1.5 in the range 10^® < Cp < 10®^ 
ergs. Two thirds of the flares considered by Crosby et al. 
(1993) fall in this range according to their Fig. 6 , so that 
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the occurrence rate of these flares can be taken as 8.7 per 
day (TVeme ~ 32.2). Eqs. imply < lO^o 

ergs and /3 f ~ 33. Assuming that the efficiency / does not 
depend on the energy, the distributions of total energy and 
of energy input per mode are identical, and N^s = 1.0. 
Eqs. ( p^ ) and (|3^) imply: 

22.4% < A < 44.6% (41) 

< cf > > 10^® ergs event“^ rnode”^, (42) 

^ ^ < 10^ modes. (43) 

< CF > 

With £^p{£) ~ Eq. (|3^) indicates that Eq. (^ 

would still be correct if the energy dependence of the effi¬ 
ciency f[£) were to select only the most energetic events. 
The size of the flare region viewed from earth is no more 
than a few arcminutes, so that all p modes with a degree 
i < ^max must receive an equal amount of energy. With 
Imax = 37 according to Wolff (1972), this corresponds to 
a set of A/” > 10^ modes, which contradicts Eq. (p^. In 
conclusion, the total acoustic energy input required from 
X-ray flares to produce the observed correlation exceeds 
their total energy. 

4-3. Coronal Mass Ejections 

4.3.1. Geometry and timing of CMEs 

Coronal mass ejections correspond to the release of 3.3 x 
10 ^^ g of matter on average, with a mean velocity of 
350 km s~^ (SMM data from 1980 and 1984-1989, Hund- 
hausen 1997). The footpoints of the outer loop are sepa¬ 
rated on average by about 45 degrees in latitude, i.e. much 
larger than an active region or flare (Harrison 1986, Hund- 
hausen 1993). This large scale structure favours low degree 
modes. According to Hundhausen (1994), the origin of the 
CME comes from a gradual build up and storage of en¬ 
ergy in a pre-ejection structure (driven by the spreading 
of magnetic field lines, the emergence of magnetic flux, or 
the shear of field lines). A finite quantity of this energy 
is then released in a catastrophic breakdown in the sta¬ 
bility or equilibrium of the stressed structure. When the 
CME is associated with an X-ray flare, the CME kinetic 
energy seems to be uncorrelated to the flare peak intensity 
(Hundhausen 1997). 

In Hundhausen (1994), the acceleration of the CME ob¬ 
served by SMM on 23rd August, 1988 was measured. Its 
launching precedes the flare event (see also Harrison 1986), 
and the acceleration to a velocity of 1000 km s“^ occurs in 
less than 10 min (see the CMEs of 17th August, 1989 and 
16th February, 1986 for similar examples in Hundhausen 
1994). 

Note that a CME was associated to the flare of 9th July, 
1996 analysed by Kosovichev & Zharkova (1998). 


4.3.2. Solar cycle variability 

While the average mass of a CME varies little from year 
to year (Hundhausen et al. 1994b), the annual variation of 
their mean velocity does not follow the solar cycle (Hund¬ 
hausen et al. 1994a). It would therefore be interesting to 
check the correlation between the high energy events no¬ 
ticed by Chaplin et al. (1995,1997) and this new indicator. 
According to Webb & Howard (1994), the occurrence rate 
of CME varies from 0.2 per day at solar minimum to 3 
per day at solar maximum. As for flares, the variation of 
the occurrence rate is enough to account for the observed 
variation of the correlation. Eq. ( p^ implies: 

0.20% < Cmin < 3.3%. (44) 

4.3.3. Energy distribution of CMEs 

Hundhausen et al. (1994a,b) showed that the distribution 
of CME velocities (10 km s~^ to 2000 km s~^) is more 
widely spread than their distribution of masses ( 10 ^"^ to 
10^® g). As a consequence, the distribution of kinetic en¬ 
ergy is spread over a much wider range than the potential 
energy. 

Hundhausen (1997) selected a sample of 249 CMEs mea¬ 
sured by SMM in 1984-1989, associated with X-ray flares. 
The spread of their distribution of kinetic energy (from 
10^® to 10®® ergs), measured from Fig. 1 of Hundhausen 
(1997) is such that /3cme ~ 20.2. Note that a smaller value 
(/3cme ~ 3.9) is obtained when considering the distribu¬ 
tion of squared velocities of 109 CMEs during the 160 
days IPHIR period in the catalogue of Burkepile & St. 
Cyr (1993). The apparent contradiction may come partly 
from the distribution of masses, but mostly from the high 
sensitivity of /3cme to the high energy events which oc¬ 
curred in 1989 at solar maximum. Moreover, fast CMEs 
are underrepresented in those statistics because the veloc¬ 
ity is measured only when the CME is visible on more than 
one coronograph image (Hundhausen et al. 1994a). As for 
flares, the distribution of acoustic energy input might be 
different from the distribution of CMEs kinetic energy, for 
example, if only a subclass of CMEs excite acoustic waves. 
In particular, some CMEs are known to be slowly accel¬ 
erated with the solar wind, while others are much faster 
than the solar wind (MacQueen & Fisher 1983). 

With /3cme = 20 and an occurrence rate of 3 CMEs per 
day {Ncme = H-l, A^eff = 0.6), IPHIR correlations in 
Eqs. ( p^ and (p^ ) imply (see Fig. ||): 

16.7% < A < 33.2%, (45) 

< Seme > > 2.4 X 10^® ergs event“® mode”®, (46) 

^ ^ < 2.8 X 10"® modes. (47) 

^ Ceme ^ 

Estimates of both A and the acoustic energy input per 
mode < Ceme > would be multiplied by a factor 2.2 if the 
value /3cme = 4 were chosen. Although the kinetic energy 
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Fig. 3. Relationship between the mean acoustic energy 
input per mode of each impulsive event and its occur¬ 
rence rate in order to produce the correlation measured 
from IPHIR or GOLF data. The kurtosis of the energy 
distribution of the source of excitation is /3ad = 20, the 
energy of the mode is 8 x 10^^ ergs and its damping time is 
Td = 3.7 days. The corresponding fractions of total power 
A = 1 and 0.3 are indicated by the continuous and long 
dashed lines. The occurrence rates of CMEs at solar min¬ 
imum and maximum are indicated by circles. 

distribution of CMEs is less accurately known than for 
flares, one can infer from the published ones (Hundhausen 
et al. 1994b, Hundhausen 1997) that the product £‘^p{£) 
is maximum at high energy. As for flares, Eq. would 
not be changed if the efficiency of acoustic wave generation 
were highest at high energy. 

Although Eqs. (p5[)-(p^ resemble Eqs. (^l|)-(^) obtained 
for flares, the following important differences should be 
noted: 

(i) the estimate for CMEs is based on their kinetic en¬ 
ergy only, which is certainly much smaller than the energy 
of the mechanism responsible for their ejection. 

(ii) the large scale of CMEs might favour the excitation 
of low degree modes, especially if their ejection mechanism 
is rooted in the convection zone {Af < 100). 

(iii) although the range of energy of CMEs (kinetic 
-I- potential ^ 8.5 x 10^° ergs on average, Hundhausen 
et al. 1994b) is the same as for big flares, their momen¬ 
tum can be one hundred times larger than the momentum 
carried by downflowing electrons produced by flares (see 
Table |^). This is favourable to a higher efficiency of the 
energy transfer to acoustic modes. 

5. Conclusion 

We have investigated the consequences of interpreting 
the observed correlations in terms of an hypothetical 


excitation mechanism, in addition to the well established 
excitation by the granules. In particular, this interpreta¬ 
tion requires impulsive events occurring no more than a 
few times per day in the IPHIR period. This has drawn 
our attention to the largest X-ray flares and CMEs at 
solar maximum. The variation of their typical energy and 
occurrence rate with the solar cycle could account for 
the variation of the correlation between the IPHIR and 
GOLF observations. 

If due to solar flares, the correlation determined from 
IPHIR data requires that at least 10“^ of the energy 
of each X-ray flares is injected into each low degree 
mode at solar maximum. Given the number of modes (at 
least 10^ due to the small scale of flares) which should 
receive the same energy as low degree modes, we are 
inclined to exclude X-ray flares regardless of the effi¬ 
ciency of the acoustic wave emission by each event. This 
reasoning, however, does not exclude the possibility of 
more energetic processes associated to flares, such as the 
restructuring of the magnetic held in the flare region men¬ 
tioned by Wolf (1972) and Kosovichev & Zharkova (1995). 

If the correlation is due to CMEs, at least 3.6 x 10“® 
of the kinetic energy of each CME should be injected into 
each low degree mode at solar maximum. This leaves two 
possibilities open: 

(i) CMEs may correlate only a few tens of low degree 
modes by injecting a few per cents of the CME kinetic 
energy into these modes. Higher degree modes {I > 10) 
would not be correlated in this case. 

(ii) CMEs could correlate more modes if significantly 
more than the observed mechanical energy of CMEs can 
be extracted from their energy reservoir and injected into 
acoustic modes. 

If all CMEs participate to acoustic emission with an 
efficiency / independent of their energy, they should be 
responsible for a fraction 16.7% < A < 33.2% of low 
degree p mode total power at solar maximum. Neverthe¬ 
less, this fraction might be significantly smaller if only 
a subset of high energy CMEs participates to acoustic 
emission. The contribution of CMEs to the low degree 
p mode power has to be reconciled with the observed 
30% decrease of their power at solar maximum. A better 
theoretical understanding of the efficiency of the energy 
transfer from a CME event to acoustic modes is needed. 

If our interpretation of IPHIR correlations is correct, 
we should be able to make the following observational 
tests: 

- detect the effect of the largest CMEs on the energy 
of low degree p modes, even at solar minimum, 

- measure the correlation Cmin > 0.2% at solar mini¬ 
mum, using longer time series and more modes than F98, 

- confirm with SOHO, during the next solar maximum, 
the correlations determined from IPHIR data, measure 
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them with a better accuracy, and determine whether 
higher degree modes are also correlated. 

Observations of the influence of CMEs on acoustic 
modes could improve our understanding of both the ex¬ 
citation of acoustic modes and the ejection mechanism of 
CMEs. 

Acknowledgements. The author thanks T. Amari, M. Tagger, 
D. Gough and D. Saundby for helpful discussions. 


Appendix A: Energy of a mode excited 
stochastically 

Let us consider a spherically symmetric, adiabatic model of 
the sun. Any perturbation of velocity v can be projected onto 
the basis of orthogonal eigenfunctions associated to the 

(supposedly real) eigenvalues tOni. The structure of the equa¬ 
tions allows us to write each component of the velocity vector 
as the product of a function of r only and a function of 8, 4> (see 
e.g. Unno et al. 1989). This function is a spherical harmonic 
for the radial velocity 


where the real dimensionless functions gf^ni^k) and h^i are de¬ 
fined as follows: 

p27T pTT 

gim{Sk) = ^ J d(l> J d(cos6')pl’"l(cos6») 


X e ’^’‘'>gk{d,dk,cl)- cl>k), 

(A.ll) 

TT 

= qim / d(l> d(cos0)P;''"'(cos6>) 


Jo Jo 


X cos{m(j>)gkid,9k,(l>), 

(A.12) 

hti = -K- [ pr^v^i{r)hk{r)dr. 

(A.13) 

Jo 



Note that the function = gf is real because of the cylin- 
dric al sym metry of the impulsion. 

Eq. (A.f) corresponds to initial conditions with zero displace¬ 
ment a.t t = tk, suitable to describe the transfer of impulsion 
due to a shock. 

Assuming that the eigenfrequencies uJni are real, we add by 
hand a damping term, with a time scale Td- After a series of 
excitations indexed by k, and neglecting the transient veloci¬ 
ties present during each event, the linearity of the equations 
allows us to write the velocity as follows: 




v„im(r) = 


«;*(r),^(r)^,n(),(r)^^ 


The eigenvectors are normalized as follows: 

pMi, 

/ = 1 , 

Jo 

= [ Pilv'nif+ 1(1+ l)\v'^i\‘^]r^dr. 


(A.l) 

Yr{e,cP). (A.2) 

(A.3) 

(A.4) 


The spherical harmonics are written in terms of Legendre as¬ 
sociated functions P;'": 


Ar,lm{t) = 2 (-l)^^^ 

X ^e cost<;„i(t - tfc). (A.14) 

tk<t 


Denoting by and the displacement and 

velocity associated to the mode the acoustic energy 

Enim of each mode is the sum of the kinetic and potential 
energies: 


1 r^* 

Enlm = 2 J + OJnl\^nlm\‘^)dMr. 


(A.15) 


Neglecting the slow damping compared to the fast oscillations 
{WniTd ~ 7 X 10^ S> Ifpr the p modes we consider), and using 
the normalization (A.3), we may write the energy as follows: 


Yr{e,cj,) = (-l)'^gi,„pl™l(cos0)e‘’"^ 

|' 2/-H (/-|m|)! ] ^ 

47r (Z + |m|)! 


(A.5) 

(A.6) 


The radial velocity perturbation Vk (r) is described by Eqs. (|^- 
(^). The angle a > 0 made by the direction d,(l> with the 
direction 6k, (j^k is defined by: 


cosa{e,ek,4’ - (l>k) = 

cosdcosOk + sind sin dk cos{(f> — (f>k). (A.7) 

By projecting the perturbation onto the basis of eigenvectors, 
we obtain: 


Anl m {t > tk) 


^nlm 


cosud„l{t - tk) 



Vfc(r) • v*(„(r)dM,.(A.8) 


2(-l) 5 *"*‘^*^ cosa;„i(t - tfc), (A.9) 

htig^m{6k)M^Vk (A.IO) 


Enlm(t) — “ ( lAniml -|- 2 

2 V 


dAnln 


dt 


(A.16) 


Using Eq. (A.14) and some algebra, we separate the two con¬ 
tributions from the waves propagating azimuthally in opposite 
directions: 


Enlmit) — E.^i.^{t) + E.^ij^{t), 
with 


/?+ — 
nlm 


e ^d a„i^e '■ 


tk<t 


JT'~ = ??"•" 

nlm nl — m' 


(A.17) 


(A.18) 


(A.19) 


Although Eqs. (A.14) and (A.16) imply E„im = Eni-m, the 
two components and E~i^, are not equal, except on av¬ 

erage. 

This separation of the components of the energy is important 
since the rotation enables us to separate these two components. 
Let us restrict ourselves to the simplest case of a solid body 
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rotation, and neg lect t he Coriolis forces. This is equivalent to 
replacing in Eq. (A. 18) the azimuthal angle </) by 0 + fit and 


4'k by 4>k + and thus we obtain Eqs. (|^-(|^. E^im is the 
energy associated to the frequency uo-ni + mfl, while is 

associated to the frequency u-ni — mfl. 

The average energy input < Cnim > of a single random excita¬ 
tion onto the mode n, I, m is proportional to the total kinetic 
energy £ of the perturbation: 


^nlm. ^ ^ ^ ■ 


(A.20) 


Appendix B: Correlation produced by a single 
excitation mechanism 

B.l. Case N 1 

In order to treat the general case, we need to establish first 
some properties of the case A <C 1. The index nlm of the 
energy is omitted in what follows, for the sake of clarity The 
mean energy < E > of the damped oscillator is directly pro¬ 
portional to the mean energy input < e > of each impulsive 


event: 





‘ 0 

"d ) dt ~ 

A 

T’ 

(B.l) 


< _£/ > ~ 

A 

^<e>, 

(B.2) 


< ~ 

A 2 

4 

(B.3) 


The energies E, E' of two modes nlm, n'l'm' satisfy: 


< EE' > 1 < ee' > 

< E X E' > ^ N < e>< e' >' 


B.2. General case 


The correlation in the general case is derived using the follow¬ 
ing theorem: if a continuous function f(x) satisfies f{2x) = 
f{x) for any value of x, and lima;^o f{x) exists, then / is con¬ 
stant. 

Let us divide the random series of excitation into two indepen¬ 
dent subset of same statistical characteristics, with an occur¬ 
rence rate A/2. Denoting by E(A/2) the energy of the mode 
excited by only one subset, one can check from Eq. that: 

< E(A) >= 2 < E(A/2) >, (B.5) 


and conclude that /(A) = E/A is a constant, which we deduce 
from the case A <C 1 (Eq. B.2): 


< E >= — <e> . 


(B.6) 


The variance of the energy can also be written as follows: 


VarE 
< E >2 


(A) 


1 Var E 
2<E>2 


{N/2) + ^. 


(B.7) 


This implies that /(A) = A(Var Ej < E — 1) is a constant, 

which we deduce from the case A ^ 1 (Eq. B.3): 


Var E _ 1 < > 

< E >2 ^ A< e >2’ 


(B.8) 


The correlation C between two modes nlm, n'l'm' is defined 
as: 

C s <EE'>-<E><E'> 

(V!irEl7(VarE')7 

Using the fact that 


< EE' > 

< Ex E' > 


(A) 


1 < EE' > 

2 < E >< E' > 


{N12) + i. 


(B.IO) 


we conclude that /(A) = A(< EE' > /(< E >< E' >) — 1) is 
a constant, which we deduce from the case A -C 1 (Eq. B.4). 
The correlation expressed by Eq. (tl2|) is then derived from 
Eqs. (p!|)-(pT|). 


Appendix C: Correlation produced by two 
excitation mechanisms 

Let us consider two independent series of impulsive events 
due to two different excitation mechanisms superimposed, in¬ 
dexed by k, characterized by the series of velocity perturbations 
and mean interpulse times The energy 

of a mode {nlm) can be decomposed into two random walks of 
Ai and A 2 steps. 


E 


.!zT 


k(l) 

“nLe ■ 


k(l) 


+ 


At) = 


E 


X2) 


- fc(2) 


fc(2) 


• (C.l) 


Let us write this equation for two modes (nlm) and {n'l'm'). 
Hereafter we simply denote the quantities specific to the mode 
{n'l'm') with a prime, for the sake of clarity. Ei and E' are 
the energies of the modes nlm and n'l'm' excited by the ex¬ 
citation mechanism {i) only {i = 1,2). The series of phases 
^k(i) ^ ^k( 2 ) ^ ^k(iy ^ ^fc( 2 )' independent. 


< E > = < El > -f < E2 >, (C. 2 ) 

Var(E) = Var(Ei) -f Var(E2) -f 2 < Ei >< E2 >, (C. 3 ) 

Var(E)C = Var(Ei)Ci-f Var(E2)C2, (C. 4 ) 


where Ci is the correlation between Ei and E(. Let us define 
the fractions A, A' of the total power of each mode contributed 
by the second mechanism as follows: 


A = 


A' = 


< E2 > 

< E > ’ 

< E2 > 

< E > ■ 


(C.5) 

(C.6) 


Let convection be the first excitation mechanism (Ci = 0), 
and let us use the index (ad) for the prop ertie s Aad, Cad o f the 
second excitation mechanism. Using Eqs. (B.f), ( p^ and (C.3) 
in (C.4), we obtain: 


C = 


< ee > 


X 1 + 


Aad < e 
A2 < e2 > 


1 + 


< e2 > 2 < e'2 > 2 

Aad <e' 

A'2 < e'2 > 


(C.7) 
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Let Pt{£) be the density of p rob ability of the additional source 
of energy as dehned in Sect. 2.5. The number of efficient exci¬ 


tations per damping time is independent of f^d- 

iVad =Nt PT(£)d£. 

Jsi 


(C.8) 


The mean energy and kurtosis of the distribution of efficient 
events are: 


£pT{£)d£ 


^ad ^ — fad 




fl^PT(£)d£ 


/3ad = 


f;yT(£)d£ £^PT(£)d£ 

£pT(£)d£y 


(C.9) 


(C.IO) 


The kurtosis is also in depe n dent of /ad. Eg. is obtained 
by incorporating Eqs. (C.8)-(C.1C) into Eq. ([2^ 
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